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ABSTRACT 

The forest of Lyman-a absorption lines seen in the spectra of distant quasars has 
become an important probe of the distribution of matter in the Universe. We use 
large, hydrodynamical simulations from the OWLS project to investigate the effect 
of feedback from galaxy formation on the probability distribution function and the 
power spectrum of the Lyman-a transmitted flux. While metal-line cooling is unim- 
portant, both galactic outflows from massive galaxies driven by active galactic nuclei 
and winds from low-mass galaxies driven by supernovae have a substantial impact 
on the flux statistics. At redshift z = 2.25, the effects on the flux statistics are of a 
similar magnitude as the statistical uncertainties of published data sets. The changes 
in the flux statistics are not due to differences in the temperature-density relation of 
the photo-ionised gas. Instead, they are caused by changes in the density distribution 
and in the fraction of hot, collisionally ionised gas. It may be possible to disentangle 
astrophysical and cosmological effects by taking advantage of the fact that they induce 
different scale and redshift dependencies. In particular, the magnitude of the feedback 
effects appears to decrease rapidly with increasing redshift. Analyses of Lyman-a for- 
est data from surveys that are currently in process, such as BOSS/SDSS-III and 
X-Shooter/VLT, must take galactic winds into account. 
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1 INTRODUCTION 

The Hi Lyman-o: forest, i.e. the multitude of Hi Lyman- 
a absorption lines seen in the spectra of distant quasars 
(QSOs), is an important cosmological observable that probes 
matter density fluctuations in the intergalactic medium 
(IGM) over a unique range of redshifts, scales and environ- 
ments. Many attempts have been made to measure physi- 
cal properties of the IGM using Lyman-a forest data. The 
two most common approaches are either based on decom- 
posing the information encoded in the transmitted flux via 
Voigt proflle fltting or treating the flux as a continuo us field 
with directly measurable sta t istical propert ies (e.g. iRauchI 
1 19981 : ICroft et al.|[l998l . |2002| : lMeiksin|[2009l ). In the second 
approach, measurement of the zero-, one-, two- or three- 
point probability distribution functions (i.e. the mean flux 
level, the flux probability distribution function, the flux 
power and bispectrum) enable a variety of physical prop- 
erties to be explored. The mean flux level for example, is 



ta-ga 

let (UV) background (e .g. iRauch et all Il997l : iTvtler et all 
I2OO4I : iBolton et all l2005h while the flux probability distri- 
bution function (PDF ) is sensitive to the thermal evolu- 

McQuin n et al.l 



tion of the IGM 



Theuns et all 



_ . 12 000 

I2OO9I : iBolton et all l2009l : iPeeples et ah 2010b a, 



The flux 

power spectrum has been used to constrain cosmological pa- 
ramete rs and the behav i our of da rk matte r (DM) at smal l 
scales jCroft et al.lll999 l: IViel et al . 2004; S eliak et al.ll2006l : 
IViel et al.ll2008l )Tthe flux bispectrum can be used to search 
for signatures of non-Gaussianities in the matter distribution 
jViel et al.ll2009l ) and wavelet techniques can also b e applied 
in order to con strain th e IGM te mperature (e.g. iMeiksinI 
l200d : iTheunTe t al. 2002i: iGarzilF et al. 201^). 

The data used for these investigations consist mainly 
of two sets of QSO spectra: the SDSS low-resolution, low 
signal-to-noise sample and UVES/VLT or HIRES/KEGK 
samples of high-resolution spectra. The characteristics of the 
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low- and high-resolution data sets are very different (the 
number of SDSS spectra is a factor ~ 200 larger than that of 
high-resolution samples, but the latter probe smaller scales 
due to the higher spectral resolution). Measurements based 
on Lyman-a forest data have reached a level of accuracy 
where an understanding of systematic uncertainties at the 
percent level or below (the magnitude of statistical errors 
associated with the SDSS sample) has become important. 

In this work we use two sets of simulations to in- 
vestigate the effects of various physical and cosmological 
processes on Lyman-a flux statistics. Firstly we use the 
suite of large cosmological hydrodynamical simulations that 
comprise the Overwhelmin gly Large Simulations project 
(OWLS; ISchave et al] I2OI0I ) to investigate the effect of 
physical processes such as metal-line cooling and galactic 
winds driven by supernovae (SNe) and active g alacti c nu- 
clei (AGN) on the flux statistics. lAltav et al] (|201ll ) and 
iTepper-Garcia et al. have shown that the OWLS 

models reproduce the observed H i column density distri- 
butions at 2; = 3 and z = Q, respectively. Secondly, we use 
sets of simulations in wh ich cosmological parameters includ- 
ing the neutrino mas s jViel et al.) |2010| ). the thermal his- 
tory of the IGM ((Viel et al.ll2009l ^ and the effect of WDM 
jViel et al.ll2"0l3 ) are varied. 

Throughout, we focus on the effects of these phys- 
ical processes on two statistics that have been shown 
to be very constraining: the flux PDF and the flux 
power spectrum. These statistics have previously been 
investigated by several authors eithe r in isolation (e.g 



McDonald et al 



2000'; 



'Jena et al."2005'; 'Becker ct al."2007l; 



Lidz et all I2OO6; Kim ct al. 2007; Bolton ct al. 200§i 
or jointly (e.g. iTheuns et all i 2000: Me iksin et all 1200 



g. lit 

IZaroubi et al.ll2006l : iDesiacgues et al. 20071 ). We will com- 
pare with the present error bars derived fro m the flux 
PDF of high-resolution UVES QSO spectra (|Kim et all 
I2OO7I) and the flux po wer measured from the SDSS data 



s 



( McDonald enilll2006l ). ^ote that ongoing surveys such as 



BOSS/SDSS-III will strongly reduce these statistical error 
bars. 

Previous simulation studies have found the effect of 
feedback on the Lyman-a forest to be relatively small, 
because the winds tend to escape into the voids, leav- 
ing the filamentary structures responsible for most of 



2OO2I: Bruscoli et all 20031: iMcDonald et al.l l2005l: 


Cen et all l2005l: iKollmeier et all 


2006: Bertone & Whitd 


20061: Borcani & Viell l2009l: 


Tornatore et all 2010l: 


Kawata & RauchI 20071: iTeDoer-Garcfa et al.l I2OI2II. How- 



ever, the precision expected from upcoming surveys implies 
that even small differences may be important. Moreover, 
previous work has generally considered only relatively 
weak feedback processes and as a result has suffered 
from "overcooling" . Inclusion of the strong feedback that 
appears to be required to reproduce observations may lead 
to much stronger effects and qualitative differences, as 
was recently found for t he low-^: matter power spectrum 
l|van Daalen et all I2OIII : ISemboloni et al.l I2OIII ). Indeed, 
the OWLS models have already been used to demonstrate 
that feedback fro m galaxy formati on changes the density 
profiles of haloes jDuffv et al.|[2O10l ) and the matter power 
spectrum (|van Daalen et al] I2OIII ) to a degree that is 
highly relevant for upcoming weak lensing experiments 



(|Semboloni et al] l201lf ). Furthermore, by using a set of 
hydrodynamic simulations, galactic feedback either in the 
form of AGN feedback and/or winds has also been shown 
to affect the properties of Damped Lyman-a s ystems and 
metal systems like CIV (e.g. iTescari et all (j2009l.l201li'l') and 
the pr operties of the intracluster medium (e.g. lFabi an et al.l 

Imi)). 

Recently, IChang et al] (|2012l ) argued that plasma insta- 
bilities cause ultra-relativisitic pairs produced by blazars to 
dissipate their energy in the IGM and that t his process dom- 
inates over photo-heating at low densities. iPuchwein et al.l 
(aoTi) demonstrated that this heating mechanism has the 
potential to significantly affect the observed Lyman-a forest 
flux statistics. Although we will confine ourselves to more 
traditional, and more localised, feedback effects, it should 
be kept in mind that other feedback mechanisms, such as 
blazar heating, may also be important. 

Here we demonstrate that the effects of outflows driven 
by feedback effects from galaxy formation may already be 
comparable to the statistical errors of past and current sur- 
veys. We find that the effect of feedback decreases rapidly 
with increasing redshift, which may make it possible to dis- 
entangle astrophysical and cosmological effects. 

The plan of the paper is as follows: in Section 2 we de- 
scribe the simulations used and the physics implemented, in 
Section 3 we compare models that include different feedback 
processes and we compare the magnitude of the feedback ef- 
fects to that of changes in the cosmological parameters and 
in the thermal history of the IGM. Finally, in Section 4 we 
present a summary and discuss future perspectives. 



2 SIMULATIONS 

In this work we use cosmological smoothed particle hydro- 
dynamics (SPH) simulations of representativ e volumes of 
the U niverse taken from the OWLS project l|Schave et al.l 
I2OIOI I. Gravitational forces and the hydrodynamic equations 
are solved with an extended version of the parallel PMTree - 
SPH code GADGET in (last described in ISpringell 120051 '). 
Our simulations assume a ACDM cosmology with parame- 
ters from the 3-yr Wilkinson Microwave Anisotropy Probe 
(WMAP) resuhfl, = 0.238, = 0.762, fib = 0.0418, h 
= 0.73, (78 = 0.74 and = 0.951. 

In addition to gravity and hydrodynamic forces, the 
simulations use sub-grid modules to follow the galaxy forma- 
tion process. The sub-grid modules used in each simulation 
include the modelling of star formation, SN feedback, the 
growth of supermassive black holes (BHs) and AGN feed- 
ba ck, radiative cooling an d chemodynamics, as described 
in ISchaye fc Dalla Vecchial ll2CI08l). iDalla Vecchia fc Schavj 



(l2008h.lBooth fc Schavd (120091 ') . I Wiersma et al.l (|2009l ') and 
|w iersma et al.l ("2009'). respectively 

Table [1] contains a list of the OWLS simulations used 
in this study. We begin here by describing our fiducial 
model (REF) before discussing the variations in the sub-grid 



^ These para meter values are largely consistent with the WMAP 
7-year results l lKomatsu et al.ll201ll) . with the main discrepancy 
being that erg is 8%, or 2.3cr, lower in the WMAP 3 data than 
allowed by the WMAP 7 data. 
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Table 1. Summary of the OWLS models used in the present analysis. The simulations that we analyse differ from one another only 
in terms of the included feedback processes, and only these parameters are noted here. Each simulation contains 512"^ particles of both 
DM and gas. A more complete description of the physics included in each simulation is present in Sec. 2. From left-to-right, the columns 
contain: (1) simulation identifier, (2) Whether the simulation treats AGN feedback, (3) Whether the simulation includes cooling through 
metal lines, (4) Whether or not the simulation contains SN winds, and for the simulations that do contain SN winds, columns (5) and 
(6) list the wind mass-loading and wind velocity. 



Simulation'-^) 


AGN(2) 


Z-cooling(^) 


SN Winds(4) 


mass-loading(^' 


wind- velocity 


REF 


X 


/ 


/ 


2 


600 


AGN 


/ 


/ 


/ 


2 


600 


NOSN.NOZCOOL 


X 


X 


X 






NOZCOOL 


X 


X 


/ 


2 


600 


WDENS 


X 


/ 


/ 






WML1V848 


X 


/ 


/ 


1 


848 


WML4V424 


X 


/ 


/ 


4 


424 


WML8V300 


X 


/ 


/ 


8 


300 



Table 2. Summary of the simulations used to probe cosmological parameters and the effect of the thermal state of the IGM. Each 
simulation is compared to a reference ACDM simulation, which is identical apart from any differences described in the fourth column. 
From left-to-right, the columns list: (1) simulation identifier, (2) Comoving box size in Mpc/h, (3) Number of baryonic particles, (4) 
Brief description of the difference between each simulation and its corresponding reference case. 



Simulation'*' 


Box Size (Mpc//i)(2) 


n(3) 

-'■''gas 


Description'*' 


WDM 


25 


512^ 


Warm dark matter IkeV vs. ACDM 


V 


60 


5123 


Neutrino mass Smi/ = 0.6 eV -I- ACDM vs. massless neutrinos -|- ACDM 


7 


20 


256^ 


T — p relation of the IGM has slope 7 = 1.6 vs. 7 = 1.3 


ris 


20 


256^ 


Spectral index ris = 0.9 vs. = 0.95 


0"8 


20 


256^ 


o-g = 0.85 vs. as = 0.8 



physics used in the other simulations. In the REF simula- 
tion, radiative cooling and heating are calculated element- 
by-element by explicitly following the 11 elements H, He, C, 
N, O, Ne, Mg, Si, S, Ca and Fe in t he presence of the Cosmic 
Microwave Background and the iHaardt k, Madaul (|200lh 
model for the UV/X-ray bac kground radiation from quasars 
and galaxies, as described in lWiersma et al] |200^Q). 

Star-forming gas is modelled by imposing a polytropic 
effective equation of state for dense particles (wh > 10~^ 
cm~^). These particles feel a pressure P oc where 
7eff = 4/3 is a polytropic index and p is the proper mass 
density of the gas. Star format ion is followed stochastically 
jSchave fc Dalla Vecchi"all2008h . with a pressure dependent 
star formation rate chosen to match the observed relation 
bet ween star formatio n surface density and gas surface den- 
sity |Kennicutt| (| 19981 ). 

The REF simulation does not include a prescription for 
BH growth and AGN feedback, but it does treat feedback 
from SNe. SN feedback is implemented by injecting ~40% of 
the energy released by Type II SNe locally as kinetic energy. 
Each newly formed star particle kicks on average = 2 of its 
neighbouring gas particles into the wind, isotropically at a 
velocity of Uw = 600 km/s. These two parameter values were 
chosen such that the peak in the global star f ormation his- 
tory a pproximately matches the observations (|Schave et al.l 
I2OIOI ). lAltav et al] l|201ll ) demonstrated that the reference 
model agrees with the observed z = 3 Hi column density 
distribution over the full 10 orders of magnitude in Hi for 
which it is measured. 

The total energy injected into SN-driven winds is pro- 
portional to r}vl,, but we have some freedom in choosing 



how the energy is distributed between the mass-loading and 
wind velocity. We therefore analyse an additional three sim- 
ulations in which the initial wind velocity and mass-loading 
are scaled such that the total energy output per unit stel- 
lar mass is the same as in REF, but the distribution across 
mass-loading and wind velocity is different. These three sim- 
ulations are termed WML1V848 (rj = 1, = 848 km/s), 
WML4V424 (tj = 4, «„ = 424 km/s), and WML8V300 
[t] = 8, w„ = 300 km/s). 

The WDENS simulation is the same as the REF sim- 
ulation but in place of using a constant wind velocity («w) 
and mass- loading (77), both parameters depend on the local 
gas density such that 

/ \ 1/8 

= eOOkms'M ^ ^'"^ A , (1) 
\0.1cm--^y 

J \ \ ^2) 

V600kms-V 

This ensures that in the densest regions winds are launched 
with the largest speeds. For gas on the star-forming equation 
of state, this scaling gives a wind speed that is proportional 
to the effective local sound speed. Note that the REF and 
WDENS models inject the same amount energy per unit 
stellar mass. The differences in results therefore result solely 
from how energy is distributed between mass loading and 
wind velocity. 

The NOSN.NOZCOOL simulation is identical to the 
REF simulation except that it does not contain feedback 
from SNe and that gas cooling rates assume zero metallic- 
ity. NOZCOOL is the same as REF except that all cooling 
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rates are calculated assuming primordial abundances. We 
can therefore isolate the effects of SN feedback by compar- 
ing simulations NOSN.NOZCOOL and NOZCOOL. 

With the exception of NOSN.NOZCOOL, all of the 
OWLS simulations analysed here contain a prescription for 
SN feedback, but the specific choices for the SN mass- 
loading and wind-velocity affect the range of halo masses 
in which SN feedback is effective. Prescriptions with high 
mass-loadings are able to eject a large amount of gas from 
objects with low star formation rates, we therefore expect 
that these prescriptions ( WML4V424 and WML8V300) will 
be able to efficiently remove gas from small galaxies at high 
redshift. These models could potentially have a large impact 
on the high-2 Lyman-a forest as it is precisely the outflows 
from these objects that have a large enough volume filling 
factor to account for observations of m etal-line absorptio n 
and hence to impact the high-z IGM fe.g. lBooth et al.]|2012l ). 

Prescriptions with high mass-loadings (and correspond- 
ingly low wind velocities) are unable to eject gas from high- 
mass galaxies because the winds are pressure-confined by 
the ambient interstellar medium l|Dalla Vecchia fc Schavel 
|2008| ) and because low wind velocities result in large ther- 
mal losses in the dense interstellar medium of massive galax- 
ies (IDalla Vecchia fc Schavell2012j ). Models that use a large 
wind velocity {WML1V848) are therefore able to eject gas 
from galaxies of higher mass. Comparison of these different 
simulations can therefore inform us how winds escaping from 
haloes of different masses affect the statistics of the IGM. 

Finally, the AGN simulation is the same as REF 
except that it includes BH growth and AGN feedback. 
These two proce s ses ar e implemented using the method of 
iBooth fc Schav3 (|2009l ) which i s, in turn, a modifie d ver- 
sion of the method pioneered bv lSpringel et all (|2005h . Seed 
BHs are inserted into the simulation by regularly running 
a halo finder and placing small BHs (rriscod = 10^ Mq) into 
every halo with mass > 10^" M0 that does not yet con- 
tain a BH. These seed BHs then grow both through merg- 
ing with other BHs and by Eddington-limited gas accretion. 
Accretion rates in low-density gas (nn < 10~^ cm~"^) are 
equal to the Bondi-Hoyle rate. In star-forming gas, this ac- 
cretion rate is boosted by a factor (nn/lO"^ cm"^)^ to com- 
pensate for the lack of a cold, interstellar gas phase and 
the finite resolution (see Booth & Schaye 2009 for a dis- 
cussion). For a given accretion rate, riiaccr, the BH growth 
rate is riiBH = (1 — er)'7iaccr, where ~ O.I is the radiative 
efficiency of the BH. 

The amount of energy available for AGN feedback is 
then given by -E = eferrfiaccrC^, where c is the speed of light 
and ef is a free parameter. Feedback is implemented by al- 
lowing the BH to build up a reservoir of accreted energy until 
it is capable of heating one of its neighbours by a temper- 
ature AT = 10* K. This energy is then injected thermally 
into one of the BH's surrounding gas particles. The param- 
eter ef = 0.15 is set by requiring that the AGN simulations 
reproduce the global BH density at z = 0. 

Each of the physics variations was run in boxes of co- 
moving size 25h~^ Mpc and 100 Mpc and contained 
512'^ particles of both gas and DM. This corresponds to an 
initial gas particle mass of 1.35 x 10^ Mq /h (8.7 x lO^M© /h) 
m the 25 (100) h-^ Mpc boxes, respectively. In Appendix 1X1 
we use a further set of simulations that use smaller box sizes 
and degraded numerical resolutions to demonstrate that our 



conclusions are robust with respect to numerical resolution 
and box size effects. 

We note that for all the OWLS models analysed in this 
paper, the "equation-of-state" (EOS) of the photo-ionised 
IGM is basically unchanged and results in the same evolu- 
tion for the parameters To, 7 that describe it. 

In order to compare the effects induced by feedback on 
the transmitted Lyman-a flux with those produced by cos- 
mological parameters or different thermal histories, we also 
consider a second set of hydrodynamical simulations. These 
simulations were also performed with the GADGET ill code 
but they include neither metal enrichment nor AGN feed- 
back and they model star formation using a simplified star 
formation criterion that converts all the gas particles with 
T < 10^ K and densities larger than one t housand times the 
mean density into stars l|Viel et al.ll2004l ). This set of sim- 
ulations is summarised in Table [2] In brief, we consider: i) 
a warm DM model (W DM) with a 1 ke V mass for the DM 
particle as presented in I Viel et al.l (|2012f ). This simulation is 
compared to a corresponding ACDM simulation (these two 
simulations were performed in a 25 comoving Mpc/h box 
with 2 X 512^ DM and gas particles); ii) two simulations that 
differ only in their choice of spectral index {us = 0.95) and 
matter power spectrum amplitude {as = 0.75), relative to 
the corresponding reference simulation (ris = 1, cts = 0.85). 
These simulations were performed in boxes of 20 comov- 
ing Mpc//i, and contain 256^ particles of both DM and 
gas; iii) two 20 comoving Mpc//i, 2 x 256^ simulations in 
which the thermal history of the IGM has different power- 
law indices for the temperature-density relation of the low- 
density photo-ioni sed gas, namely T = To(l -I- 5)^~'~^\ with 
7 = 1.3, 1.6 as in IViel et al.l (2009); iv) a simulation of 60 
comoving Mpc/h with 512'^ massive neutrinos of total mass 
Em^ — 0.6 eV and 512^ DM particles compar ed with a 
simulation in which the neutrinos are massless l|Viel et al.l 
I2OIO). 

The mass per gas particle is 1.5 x 10*^ Mq for the WDM 
simulation, 5.8 x 10® Mq for the neutrino simulation, and 
6.1 x 10® Mq for aU the other cases. This second set of simu- 
lations is not homogeneous, so we will only use it to highlight 
the relative differences with respect to the matched reference 
simulation that corresponds to each simulation (erg = 0.85, 
Us = 0.95, 7 = 1.3 and Em^ = 0). The resolution of each of 
these simulations is sufficient to obtain converged results for 
the Lyma n-a PDF at z < 3 (|Bolton et al.ll2008l : IViel et all 
12009. : iCalura et al.ll20ia ). 



3 FLUX STATISTICS 

In this Section we focus on two different Lyman-a flux statis- 
tics that are widely used for cosmological purposes and for 
constraining the physical state of the IGM: the flux PDF 
and the flux power spectrum. From each simulation we ex- 
tract a set of 1000 random lines-of-sight (LOSs) at z = 2.25, 
z = 3 and z = 4 and compute the transmitted flux along 
each of the LOS in redshift space. The synthetic spectra were 
convolved with a Gaussian with FWHM 6.6kms~^ and re- 
binned onto pixels of size 5kms~^. We generally report our 
results in terms of ratios or relative differences with respect 
to the corresponding reference model at the same resolution, 
redshift and box size. 
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Figure 1. Flux PDF (left) and flux power (right) for two different cosmological boxes (100,512) and (25,512) represented by the dotted 
and solid lines, respectively. We plot the flux PDF and flux power for the AGN (red curves) and REF (black curves) simulations at 
z = 2.25. 
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Figure 2. Effect of AGN feedback on the flux PDF (left) and flux power (right) for two different cosmological boxes (100,512) and 
(25,512) represented by the black, dotted and red, solid lines, respectively. We plot the fractional differences between the flux PDF and 
flux power of the AGN and REF s imulations at 2 = 2.25. The blue, s olid lines bracket the la statistical observed errors from UVES/VLT 
at 2 = 2.07 PDF l|Kim et al.ll2007l ') and SDSS flux power at 2 = 2.2 l|McDonald et al.ll200d l. 



In all the plots shown, we have scaled the optical depth 
of each simulation to reproduce the same mean flux level as 
meas ured from o bserved UVES/VLT high resolution spec- 
tra of iKim et af] (|2003), namely: t^b = 0.0023 x (1 + zf-^^ 
and TeS= —In < F >. This procedure is justified because 
the intensity of the ionising background radiation is poorly 
constrained. Note, however, that the original mean fluxes 
differ usually by less than one percent between the different 



simulations at all redshifts. In terms of overall scaling factors 
for the optical depths, the simulations differ by less than 7% 
(the scaling factors are also very close to one except for the 
z — 4 outputs). 

We do not consider the impact of resolution and noise 
on our mock quasar spectra because these factors do not 
alter any of our concl usions with regards to the flux PDF or 
power spectrum (e.g. iBolton et al.ll2008l ). 
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Figure 3. Effect of metal cooling, SN feedback and galactic winds in terms of flux PDF (left) and flux power (right) for the (25,512) 
simulations at z = 2.25. We plot the fractional differences between the flux PDF and flux power of the REF simulation and the following 
ones: without SN feedback and with primordial cooling (green, solid), with strong SN feedback (black, dotted), with SN feedback and 
with primordial coolin g (orange, dashed). The blue, soli d lines bracket the Icr statistical observed errors from UVES/VLT at z = 2.07 
PDF l lKim et al.ll2007^ and SDSS flux power at z = 2.2 jMcDonald et al.ll2006l) . 




Figure 4. Effect of different galactic wind prescriptions in terms of flux PDF (left) and flux power (right) for the (25,512) simulations 
at z = 2.25. We plot the fractional differences between the flux PDF and flux power of the REF simulation and the following ones: 
WML8V300 (green, solid), WML1V848 (black, dotted), W ML4 V424 (oran ge, dashed curve). The blue, sol id lines bracket the l a 
statistical observed errors from UVES/VLT at z = 2.07 PDF l lKim et al.ll2007^ and SDSS flux power at z = 2.2 llMcDonald et al.ll20od) . 



3.1 Redshift z = 2.25 

3.1.1 Comparison of OWLS models 

In Figures [T] and [5] we show the effect of AGN feedback on 
the z = 2.25 flux PDF (left panel) and flux power (right 
panel) for two different box sizes (100,512) and (25,512) 



represented by dotted and solid lines, respectively. Figure 
[T] shows the actual PDF and flux power in the two cases, 
while in Figure[2]we plot the fractional differences. The blue 
lines bracket the la error b ars obtained from the flux PDF 
measurement at z = 2.07 of iKim et al.l (l2007f) and the flux 
power spectrum measurement at z = 2.2 of McDonald et al.l 
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Figure 5. Effect of cosmological parameters and thermal state of the IGM on the z = 2.25 flux PDF (left) and flux power (right) 
compared to the AGN effect (fractional differences are shown). We show the effect of having: a lower ag value (erg = 0.75) compared to 
the reference case erg = 0.85 (black dash-dotted); a lower value (nB=0.9) compared to the reference ns=0.95 value (blue, dotted); a 
higher 7 value (7 = 1.6) compared to the reference 7 = 1.3 value (green, solid); a thermal WDM relic of 1 keV compared to the ACDM 
case (orange, dash-triple-dotted); a neutrino simulation with Sm,i, = 0.6 eV (cyan, dashed) and including AGN feedback (red, thick 
solid). See the text for a more complete description of the simulations used. 



For the flux PDF the agreement between the two 
simulations is at the 2-3% level. AGN feedback results in an 
increase of regions with high transmissivity, F > 0.5, at the 
expense of low-transmissivity regions. This is due to the fact 
that, although the simulations share the same T — p density 
relation for the photo-ionised, low-density IGM, the gas is 
distributed differently in the T—p plane in the two scenarios, 
with AGN feedback increasing the amount of low-density gas 
with T > 10^ K relative to the REF simulation. In this case 
the effect induced is comparable to the Icr (statistical) error 
bars of the measured flux PDF. In terms of flux power, we 
find a scale-dependent suppression in the AGN case which 
can reach the 10% level at the largest scales probed and is 
~ 5% at fc = 0.6 /i/Mpc. We find that AGN feedback is likely 
to suppress the flux power at 2; « 2 by a scale-dependent 
amount which is slightly larger than the Icr statistical error 
bars at fe = 0.2 - 0.6 fe/Mpc . 

It is important to note that the flux power spectra of 
the two AGN simulations do not match perfectly due to 
the large difference in the numerical resolution of the two 
runs (a factor 64 in the particle mass). Our implementa- 
tion of BH growth is sensitive to the numerical resolution 
for BHs with masses similar or smaller than the gas parti- 
cle mass l|Booth fc Schavell2009l ). While the low-resolution 
AGN model reproduces optical an d X-ray observatio ns of 
low-redshift groups of galaxies (Mc Carthy et al.]|2010l ). the 
same may not be true for the higher-resolution run, which 
was too expensive computationally to run past z = 2. 

In Figure [3] we show three different effects on Lyman- 
a flux statistics: the absence of metal cooling (solid green), 
the presence of more efficient SN-driven galactic winds (dot- 
ted black) and the absence of metal cooling and SN feedback 



(dashed orange). This Figure refers to z = 2.25. Turning 
off metal-line cooling has very little effect, but if in ad- 
dition SN feedback is turned off, then there are substan- 
tial differences. Turning off SN feedback suppresses high- 
transmissivity [F — 0.6 — 0.8) regions, so the effect on the 
flux PDF is consistent with that of turning off AGN. Increas- 
ing the strength of SN-driven winds in dense gas ( WDENS) 
has a similar effect on the ffux PDF as including AGN feed- 
back. The effects of feedback processes on the flux power are 
more confusing: turning on SN feedback (i.e. moving from 
orange, dashed to green, solid) enhances the power on large 
(k < 2 /i/Mpc ) and suppresses it on small scales, which 
goes in the opposite direction as the effect of turning on 
AGN (see Fig. [2]). Enhancing the strength of winds in dense 
gas has very little effect on the flux power except on very 
small scales. 



In Figure |3] we plot the fractional differences 
between REF and the three different wind models 
WML1V848,WML4V424 and WML8V300 in terms of the 
flux PDF (left) and flux power (right) at ^ = 2.25. For the 
flux PDF the trends are consistent with the ones shown 
in Figures [5] and |31 increasing the efficiency of feedback in 
higher-mass galaxies (i.e. using higher wind velocities) in- 
creases the fraction of high-transmissivity regions at the ex- 
pense of the low-transmissivity regions. For the flux power 
the trend is consistent with that of Fig. [S] more efflcient 
winds decrease the power on small scales, but enhance it on 
the larger scales. As noted above, this is qualitatively differ- 
ent from the effect of AGN feedback, which suppresses the 
power on large scales (see [2]). 
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3.1.2 Other simulations 

In Figure[5]we compare the effect of AGN feedback to effects 
induced by clianging some cosmological and astropliysical 
parameters. In this latter case, we use pubhshed results that 
have been obtained from other sets of simulations. However, 
even for these cosmological and astrophysical parameters 
we compare with a reference case which is reasonably close 
to the OWLS REF simulation. We can see that at z « 2 
the effect of AGN feedback on the flux PDF is comparable 
to reducing the value of ag, by 0.1 and that is larger than 
decreasing by 0.05. Changing 7, the power-law index of 
the temperature-density relation T = To(l + 5y~^ , from 
1.3 to 1.6, or using a WDM matter power spectrum with 
thermal mwDM = 1 keV particles, produces effects that are 
larger than that of the feedback and for the case of WDM 
they are also of opposite sign. 

It is also interesting to compare these models in terms 
of the flux power (right panel): the effect of including AGN 
feedback is quite scale dependent, at least at scales k < 
1/i/Mpc, while the other effects could be modelled as an 
overall suppression or increase of power. The magnitude of 
the AGN feedback effect at the largest scales is comparable 
to that induced by a different thermal history and that pro- 
duced by a WDM relic. The neutrino induced suppression 
of the flux power is of the order of 7-8% on the largest scales 
(note, however, that the neutrino suppression is computed 
in a model in which the r.m.s. amplitude of the matter power 
is normalised at the cosmic microwave background scale). 



3.1.3 Keeping the temperature- density relation fixed 

The differences in the flux statistics between simulations us- 
ing different feedback prescriptions could be the result of 
changes in the distribution of the gas density, changes in 
the ratio of photo-ionised and coUisionally ionised gas (i.e. 
the fraction of the gas with T ^ 10^ K), or it could be due 
to changes in the T — p relation of the photo-ionised gas. As 
the latter relation ca n be measured from the widths of the 
absorption lines fe.g. ISchave et aDll999l l. it is important to 
check the role of changes in the thermal state of the diffuse, 
photo-ionised gas. Below we will follow convention and refer 
to the T — p relation of the photo-ionised as the "equation 
of state" (EOS), even though the true equation of state is 
that of an ideal gas. 

We re-analysed a number of simulations after imposing 
the same EOS in a number of different ways. In each case, 
the EOS is imposed a-posteriori by changing the tempera- 
tures of (a subset) of particles and the ionization fractions 
and spectra are then recalculated. As before, the optical 
depths are re-scaled so th at the mean flux d ecrement agrees 
with the observations of iKim et all (|2007l ) for all models. 
We consider four different cases: 

• Full simulation: the thermal state of the gas predicted 
by the simulation is used, i.e. no EOS is imposed (red curve). 

• Single power-law: an EOS of the form T = ro(l + 
^^(7-1) jg imposed on all gas at all densities, with parameters 
To = 4.66 and 7 = 1.4 (green, dashed curve). This EOS fits 
the p — T relation in the low-density IGM, but overpredicts 
the temperatures at high density, where radiative cooling 
is more important than the adiabatic cooling due to the 



Hubble expansion. We include this case because a single 
power-law is a common theoretical assumption. 

• Double power-law/low-temp only: gas at high tem- 
peratures (r > 10'' K; this gas is mostly coUisionally ionised) 
retains its original thermodynamic state. At low temper- 
atures, all gas is placed onto an EOS. At low densities 
{5 < 20) the EOS is the same as described above. At high 
densities {S > 20), the EOS has a slope of 7 = 0.8. This 
double power-law EOS provides a good fit to the T — p re- 
lation of the photo-ionised gas in the simulations (orange 
continuous curve). 

• Double power-law/no hot gas: this model is the 
same as the one described directly above, but the density 
of the hot gas (T > 10^ K) is set to zero so that it does 
not contribute to the absorption (black dotted curve). This 
allows us to check the effect of absorption by hot gas. 

In Figure|n]we show the fractional differences in the flux 
PDF and the flux power between the z = 2.25 outputs of 
the AGN and REF runs, after both simulations have each 
of the EOSs described above imposed. Note that the red, 
solid curve ("full simulation") is identical to the one shown 
in Figure [H 

We can note the following trends. If we impose a "sin- 
gle power- law" EOS, then the difference between the AGN 
and REF models is completely different than if we use the 
predicted temperatures ("full simulation"). This is because 
a single power-law strongly overpredicts the temperature of 
high-density gas, which now becomes coUisionally ionised 
and hence mostly invisible in HI. Moreover, the hot, coUi- 
sionally ionised low-density gas is now given a low temper- 
ature, which increases its HI fraction. Clearly, imposing a 
single power-law EOS on all gas results in large changes to 
the flux statistics. 

The agreement with the "full simulation" is, however, 
very good if we use the "double power-law/low-temp only" 
EOS that results in colder temperature for the high-density 
gas and which leaves the gas with T > 10^ K unchanged. 
Imposing this EOS hardly changes the difference between 
the AGN and REF models. The situation is very different 
if we impose the double power-law EOS also on gas that 
originally had a temperature T > 10''' K (not shown). In 
that case the result is more similar to the "single power- 
law" model than to the "double power-law/low-temp only" 
model. 

We can draw two important conclusions from these 
tests. First, the differences in the flux statistics between the 
AGN and REF simulations is not due to changes in the EOS 
of the photo-ionised gas, because the difference between the 
two models is nearly identical if we impose the same, realis- 
tic EOS on both simulations. Second, the difference between 
models AGN and REF is not only caused by changes in the 
density distribution, because we obtain completely different 
results if we impose the double power-law EOS also on gas 
with T > 10^ K. The fact that AGN feedback heats gas to 
temperatures sufficiently high for the gas to become highly 
coUisionally ionised is very important. 

Finally, we note that the gas with T > 10^ K is not 
completely invisible in HI. We can see this by comparing 
the cases "double power-law; low-temp only" and "double 
power-law/no hot gas". The differences are always below 
the 4% level in terms of both flux power and flux PDF and 
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Figure 6. Comparison between the AGN and REF runs at 2 = 2.25 after imposing an identical T — p relation (EOS) on the gas in the 
two models. The red thick curve refers to the standard case in which no temperature rescaling is performed ("full simulation"); "single 
power-law" refers to a case in which all the gas is placed on the EOS that is appropriate for the low-density, photo-ionised IGM (green, 
dashed curve); "double power-law/low-temp only" refers to a case in which gas with T > 10® is left unchanged and lower temperature 
gas is placed on a double power-law EOS with separate slopes for gas with densities above and below 10"'* cm~^ {& = 20) (orange, 
continuous curve); "double power-law/no hot gas" refers to the same case as "double power-law" except that the HI fraction of gas with 
T > 10® K is set to zero (black dotted curve). 



of the order of 1-2% at the wavenumbers considered in the 
SDSS flux power analysis. Thus, excluding the hot gas at 
temperatures above 10® K has a small but non-negligible 
impact on the flux statistics. 

3.2 Evolution 

Up till now we have focused on redshift z — 2.25, but the 
flux statistics have been measured over a substantial red- 
shift range. At prese nt, PDF measurem ents span the redshift 
range z = 2 — 6 (e.g. iBecker et al.ll2007l ) an d flux power mea- 
surements cover the range z = 2.1 — 4.4 (|Croft et al.ll2002l : 
IViel et al.ir2004l : iMcDonald et al.ll2006l ). The large redshift 
range may allow one to lift the degeneracies between cos- 
mological and astrophysical parameters and results in tight 
constraints on the recovered parameters. In this section we 
will therefore investigate the evolution of the impact of the 
physical effects by comparing predictions for z = 3, 4 with 
the results for z — 2.25 that we already showed in the pre- 
vious sections. 

We caution that the predictions for 2 = 3 and especially 
for z = 4 are probably not fully converged with respect 
to the numerical resolution, except for the flux power on 
large scales (see Appendix As we will limit ourselves to 
comparing different physical models at the same resolution, 
even predictions that are not fully converged are informa- 
tive. However, we must keep in mind that the comparisons 
may be affected by the lack of convergence. In particular, 
the impact of winds driven by low-mass galaxies may well 
be underestimated when the star formation rate in these 
objects is limited by the finite resolution of the simulation. 

In Figure [7] we compare the AGN and the REF sim- 
ulations at z = 2.25, 3, and 4. The induced suppression of 
the flux power is 3% smaller at z = 3 than at 2 = 2.25, 



while at 2 = 4 the models agree since AGN activity, which 
is restricted to relatively massive galaxies in our AGN simu- 
lation, is rare at such high redshifts. At high redshift there is 
a substantial difference in the flux PDF for very high trans- 
missivity regions {F > 0.85), but this would be difficult to 
observe due to continuum fitting errors. 

We investigate the redshift dependence of the impact 
of SN feedback, for the subgrid recipe employed by the 
REF simulation, by comparing the NOZCOOL and the 
NOSN.NOZGOOL models at 2 = 2.25, 3, and 4 in Figure H 
On large scales, k < 2ft/Mpc, the effect on the flux power is 
nearly identical at 2 = 2.25 and 3 and at 2 = 4 the difference 
between the two models is less than 2%. On smaller scales 
the evolution is stronger. Interestingly, by 2 = 4 the effect of 
SN feedback only exceeds 2% for k > 10/i/Mpc. The impact 
of SN feedback on the ffux PDF actually increases with red- 
shift for very high transmissivity regions {F > 0.85), but as 
noted before, continuum fitting errors, which also increase 
with redshift, in any case limit the usefulness of such regions. 

Figure [5] shows the evolution of the difference between 
simulations WML8V300 and REF. This comparison is in- 
teresting, because for low-mass galaxies the winds in model 
WML8V300 are more efficient than for REF. This is op- 
posite to AGN, because AGN feedback only drives strong 
outflows in massive galaxies. Nevertheless, by 2 = 4 the dif- 
ference between WML8V300 and REF is below 2% for both 
the flux PDF and the flux power, except again for the fre- 
quency of very high transmissivity regions and the power 
on small scales. The WDENS simulation (not shown) has 
instead very little dependence on redshift in terms of ffux 
power. 

We conclude that the three different feedback scenar- 
ios have distinct signatures in the {k, 2) parameter space 
and that the feedback effects become weaker with increas- 
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Figure 7. Comparison of tlie effect of AGN feedback on the flux PDF (left) and flux power (rigtit) at redshifts z = 2.25 (solid, red), 
z = 3 (dotted, black), and z = 4 (daslied, blue). Except for the frequency of very high transmissivity regions, the effect of AGN feedback 
decreases with redshift. 
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Figure 8. Comparison of the effect of SN feedback, using the subgrid recipe employed by the REF model, on the flux PDF (left) and 
flux power (right) at redshifts z = 2.25 (solid, red), 2 = 3 (dotted, black), and 2 = 4 (dashed, blue). Except for the frequency of very 
high transmissivity regions, the effect of SN feedback decreases with redshift. 



ing redshift, except for the flux PDF in high transmissivity 
regions. Although many of the feedback effects are too weak 
to be constrained by present data sets, it will be important 
to consider them properly when the statistical error bars are 
reduced further, as will be the case for ongoing surveys such 
as BOSS/SDSS-III. 

4 CONCLUSIONS AND DISCUSSION 

In this paper we have investigated the effect of feedback from 
galaxy formation on two widely used Lyman-a flux statis- 
tics: the PDF and the power spectrum of the normalised 
flux. We have used state-of-the-art hydrodynamical simu- 
lations from the OWLS project that incorporate most of 



the physical processes that may affect the properties of the 
IGM. The suite of OWLS runs, of which we have used only 
a subset, is particularly well suited for our purposes as the 
OWLS project was designed to test the effect of uncertain 
(subgrid) processes by varying one parameter at a time in 
simulations that start from identical initial conditions. 

The impact of feedback on Lyman-a statistics is gener- 
ically small because galactic winds tend to escape into 
the voids, leaving the fllamen ts responsible for th e Lyman- 
Q absorption largely in tac t (|Theuns et al.l [20021 ). Indeed, 
iTepper-Garci'a et al.l (|2012l ) demonstrated this to be true 
for the OWLS models by comparing the statistics of Voigt 
profile decompositions of synthetic spectra at low red- 
shift. However, even small differences may be important as 
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Figure 9. Comparison of tlie difference in the flux PDF (left) and flux power (rigiit) predicted by simulation WML8V300 and REF at 
rcdshifts z = 2.25 (solid, red), 2 = 3 (dotted, black), and z = A (dashed, blue). Compared with REF, the prescription for SN feedback 
employed by model WML8V300 results in stronger outflows in low- mass galaxies. Except for the frequency of very high transmissivity 
regions, the difference between the flux statistics predicted by the two simulations decreases with redshift. 



the statistical errors on published data, such as the flux 
PDF me asured from high-resolution UVES QSO spectra 
(|Kim et a l. 2007) and the flux power measured from the 
SDSS data (|McDonald et all f20OT ). are already as small 
as 5%. Moreover, upcoming observational campaigns such 
as the BOSS/SDSS-III Lyman-a survey (jEisenstein et al] 
I2OIII : ISlosar et al.ll201ll ) and the X-S hootcr/VLT medium 
resolution spectrograph (jVernet et al.. 2011) are expected to 
reduce these error bars substantially. 

We investigated the effect of metal-line cooling, a num- 
ber of different implementations of outflows driven by feed- 
back from star formation, and AGN feedback by comparing 
different OWLS models to the OWLS reference simulation 
(REF). For each redshift the spectra drawn from each sim- 
ulation were re-scaled to the same, observed mean flux, thus 
removing a potential source of differences between the mod- 
els (this procedure is justified because the intensity of the 
ionising background radiation is poorly constrained). We fo- 
cused on redshift z = 2.25 but also studied the evolution 
from z = 4. 

Our main conclusions can be summarised as follows: 

• Metal-line cooling has a much smaller effect on the flux 
statistics than galactic outflows. 

• AGN feedback suppresses the flux power on large scales 
(fc < 1 /i/Mpc ). The effect increases in strength with de- 
creasing wave number. At z — 2.25 it reaches ~ 10% on the 
largest scales we can measure (fc ~ 0.1 /i/Mpc ). The flux 
PDF is changed at the 5% level at this redshift. 

• When using the implementation of REF, winds driven 
by SNe have an effect that is qualitatively opposite to that of 
AGN feedback. The changes in the PDF and the flux power 
spectrum caused by SN feedback are slightly smaller than 
for AGN feedback. 

• More efficient, but energetically feasible, winds from 
low-mass galaxies have a strong impact on the fiux statis- 
tics. At 2 = 2.25 the effect is similar, both qualitatively and 



quantitatively, as that of AGN feedback. At higher redshift 
the impact of such winds changes and is larger than that of 
AGN feedback. 

• Overall, the effects induced by feedback from AGN and 
star formation are of the order of the present la statistical 
error bars, comparable to those induced by cosmological and 
other astrophysical parameters. 

• All the feedback mechanisms studied here have a sub- 
stantially smaller impact on the Lyman-a forest at higher 
redshift. For winds from low-mass galaxies we can, however, 
not rule out that the evolution would be weaker for higher 
resolution simulations. 

• The feedback effects are not due to changes in the 
temperature-density relation of the diffuse, photo-ionised 
gas. Instead, they are due to changes in the density dis- 
tribution of the gas and to changes in the fraction of hot, 
coUisionally ionised gas at a fixed density. 

We have shown that the effects of feedback from galaxy 
formation on Lyman-a flux statistics are not negligible. In 
particular, ai z — 2.25 they are comparable to present ob- 
served statistical uncertainties. As upcoming surveys will 
strongly reduce these statistical errors, it will become neces- 
sary to properly account for these effects in future analyses 
of Lyman-a forest data sets. Although winds driven by feed- 
back from AGN and star formation are difficult to model 
from first principles, progress can be made by using the fact 
that their effect on the Lyman-a forest will generically have 
a different scale and redshift dependence than changes in 
the cosmological model. 
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We believe that numerical errors induced by poor resolution 
are likely to be smaller than the values quoted above when 
the (25,512) simulations are considered (we would require a 
higher resolution simulation to confirm this). The fact that 
we are interested in relative differences between different 
models using the same resolution will largely mitigate any 
resolution eflFects. An exception may be the importance of 
feedback from low-mass galaxies, which we will underesti- 
mate if the star formation rate in low-mass haloes is limited 
by the finite resolution. This may be the case at higher red- 
shifts. 



APPENDIX A: BOX-SIZE AND RESOLUTION 
EFFECTS 

In this appendix, we consider box-size and resolution eflects 
for the reference (REF) runs. In Figure ET] we quantify the 
effects that a different box size has on the flux PDF (left 
panel) and on the flux power (right panel) by considering 
three different boxes of linear comoving size 100, 50 and 
25 Mpc at fixed resolution and at 2 = 2.25. From the 
Figure we conclude that the 50 and 100 Mpc boxes agree 
well for k ^ 1/i/Mpc. The fact that large scale modes are 
affected is not surprising, but this issue is not critical for us 
since we are interested in the differences of various physical 
models relative to the reference model run m the same box 
size from the same initial conditions. 

Figure rX2l auantifies resolution effects by comparing two 
reference runs: the (25,256), represented by the continuous 
red curve, and the (25,128), represented by the dot-dashed 
red curve, to the reference run at higher resolution (25,512). 
Note that the mass (spatial) resolution varies by factors of 
8 (2) between each pair of simulations. At z = 2.25 the flux 
PDF is converged to better than 5% over the whole flux 
range and the flux power is also converged at the 3% level 
up to fc = 3/i/Mpc between the (25,256) and (25,512) runs. 
At higher redshifts, 2 = 3, 4 (dashed, orange and dotted, 
green curves, respectively) we compare the (25,256) and the 
(25,512) simulations and find that the convergence is worse. 
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Figure Al. Effect of different box size at fixed resolution for the REF simulations a,t z = 2.25. We plot the difference of the (25,128) 
(black, dotted) and (50,256) (green, solid) boxes relative to the (100,512) simulation, respectively (the flux power results are only shown 
for the range of wavenumbers in which there is overlap between the different boxes). 




Figure A2. Effect of numerical resolution for the reference (REF) simulations at 2 = 2.25, 3, 4. We plot the relative difference of the 
(25,256) and the (25,512) simulation as the red solid (z = 2.25), dashed orange (z = 3) and green dotted curves [z = 4). We also show 
the ratio of the (25,128) and (25,512) boxes at z = 2.25 as a red dot-dashed curve. 
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